Desenvolvimento de modelo de primeiros princípios para simulação dinâmica de uma coluna de destilação contínua

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Química e Biológica

Simulação Dinâmica do Centro de Dia e Apoio à Terceira Idade de Leça do Balio

Pela importância que os edifícios têm na utilização de energia, a avaliação do seu desempenho energético é de grande relevância, uma vez que, em grande parte, passa por estes a concretização das metas europeias definidas para 2020, no que concerne à diminuição da utilização de energia. Tendo em conta que os edifícios representam 40% do consumo de energia total, e estando o sector em expansão, esta realidade obriga a uma procura de soluções integradas de arquitetura e engenharia que promovam a sustentabilidade dos edifícios. Foi efetuado um estudo num edifício constituído por dois corpos, um mais antigo que funciona como centro de dia e um mais recente que funciona como lar, localizados no concelho de Matosinhos, onde se identificaram os pontos de maior consumo energético, para os quais foram sugeridas alterações no sentido de baixar os custos com a factura energética. Nesta dissertação foi utilizado um software de simulação dinâmica para avaliação do comportamento térmico do edifício nas condições atuais e, posteriormente, foram simulados outros cenários com alterações ao nível da envolvente térmica dos edifícios e dos seus sistemas técnicos, que permitiram identificar algumas medidas de melhoria de eficiência energética. As medidas de melhoria sugeridas implicam uma redução energética, ao nível do consumo de água quente sanitária, consumo de gás natural e electricidade. De entre essas medidas, e com um payback inferior a 8 anos e meio, destacam-se a instalação de redutores de caudal, a substituição da caldeira e da bomba de recirculação, a instalação de painéis solares térmicos e a redução da quantidade de lâmpadas.

Análise de sensibilidade na simulação dinâmica de um grande edifício escolar existente

Dado o panorama de conservação de energia a nível nacional e mundial, torna-se hoje em dia muito importante, que seja possível controlar e estimar o consumo energético nos edifícios. Assim, atendendo à actual problemática energética e ao crescente consumo energético nos edifícios, é importante parametrizar, avaliar e comparar este consumo. Neste sentido, nas últimas décadas, têm sido efectuados desenvolvimentos técnicos, quer ao nível do equipamento de campo para efectuar monitorização e medição, quer ao nível da simulação dinâmica de edifícios. Com esta dissertação de mestrado, pretendeu-se efectuar a simulação dinâmica de um edifício escolar existente a funcionar em pleno, e efectuar uma análise de sensibilidade relativamente ao grau de variação dos resultados obtidos através da aplicação de dois programas de cálculo térmico e energético. Foram utilizados, o programa VE-Pro da IES (Integrated Environmental Solutions) e o programa Trace 700 da TRANE. Ambos os programas foram parametrizados com os mesmos dados de entrada, tendo em atenção as opções de simulação disponibilizadas por ambos. Posteriormente, utilizaram-se os dados retirados da simulação para calcular a classificação energética no âmbito do sistema de certificação energética (SCE), através de uma folha de cálculo desenvolvida para o efeito. Foram ainda consideradas várias soluções de eficiência energética para o edifício, com vista a poupanças reais de energia, tendo sempre atenção ao conforto térmico dos ocupantes. Dessas medidas fazem parte, medidas relacionadas com a iluminação, como a substituição da iluminação existente por luminárias do tipo LED (Light Emitting Diode), soluções de energias renováveis, como a instalação de colectores solares para aquecimento das águas quentes sanitárias, e painéis fotovoltaicos para produção de energia, bem como medidas ligadas aos equipamentos de climatização. Posteriormente, recalculou-se a classificação energética afectada das melhorias. Os resultados obtidos nas duas simulações foram analisados sob o ponto de vista do aquecimento, arrefecimento, ventilação, iluminação e equipamentos eléctricos. A comparação das duas simulações para cada parâmetro acima referido, apresentaram variações inferiores a 5%. O desvio maior verificou-se na ventilação, com o valor de aproximadamente 4,9%. Transpondo estes resultados para o cálculo do IEE (Índice de Eficiência Energética), verificou-se um desvio inferior a 2%.

Estimação e simulação dinâmica de um modelo econométrico para o mercado de alumínio no Brasil

O presente estudo pretende analisar o mercado brasileiro de alumínio primário, através da elaboração e simulação dinâmica de um modelo econométrico, com capacidade de explicar satisfatoriamente as relações típicas das variáveis econômicas consideradas neste mercado, conforme são postuladas pela teoria da empresa em competição imperfeita.

Especificação e simulação dinâmica de um modelo macroeconômico de pequeno porte para a economia brasileira

None

Análise do setor da agências de internet no Brasil : uma aplicação de simulação dinâmica de sistemas

Esta dissertação tem por objetivo analisar a participação de mercado dos diferentes grupos estratégicos que compõem o setor de agências de internet no Brasil, bem como as mudanças decorrentes de alterações nas variáveis críticas de concorrência, tais como preço, investimento em propaganda, investimento em tecnologia e disponibilidade de mão-de-obra.

Desenvolvimento de framework para análise e simulação dinâmica de sistemas elétricos de potência

Esta tese apresenta uma metodologia flexível orientada a objetos (OO) para a aplicação no projeto e implementação de sistemas de software utilizados na realização de estudos dinâmicos de sistemas elétricos de grande porte. A metodologia OO proposta objetiva tornar mais simples o desenvolvimento, a atualização e a manutenção de complexos sistemas de software para estudos de transitórios eletromecânicos em sistemas elétricos de potência. Os requisitos de usuário são mapeados para um conjunto de classes básicas, as quais são usadas para efetuar a modelagem de dispositivos dinâmicos tais como geradores elétricos. Para avaliação da metodologia foram realizados dois estudos de casos. No primeiro estudo caso o Framework foi aplicado na simulação das unidades geradoras da Usina Hidrelétrica de Tucuruí. Os resultados da simulação foram comparados com medições obtidos em ensaios no campo e mostrou a boa performance do Framework na reprodução dos fenômenos eletromecânicos desta usina de grande porte. No segundo estudo de caso, por outro lado, o Framework foi aplicado na modelagem de um sistema de geração fotovoltaico (PV) com seu sistema de Rastreamento da Potência Máxima (MPPT). O controle MPPT foi implementado usando técnicas digitais. Os resultados das simulações demonstram a performance do Framework na modelagem do sistema de controle de corrente, assim como no controle MPPT, dos sistemas de geração PV.

Modelagem e simulação da dinâmica do arrastamento por ar

Os compostos orgânicos voláteis constituem uma fonte vulgar de contaminação da água subterrânea, a qual pode ser eliminada pela tecnologia do arrastamento por ar (air stripping) em colunas com enchimento desordenado e utilizando fluxos das fases em contra-corrente. Propõe-se neste trabalho uma nova metodologia de dimensionamento destas colunas, para qualquer tipo de enchimento e de contaminante, onde não há necessidade de se arbitrar nenhum diâmetro, onde se evita o recurso a ábacos experimentais e onde o regime hidráulico conveniente é seleccionado à partida. O procedimento proposto foi algoritmizado e convertido num programa em linguagem C++. Para verificar e testar não só o dimensionamento mas também o comportamento teórico estacionário e dinâmico construiu-se de raiz uma coluna experimental. Seleccionou-se como contaminante uma solução de clorofórmio em água destilada. A experimentação permite, ainda, corrigir o coeficiente de transferência de massa global teórico estimado pelas correlações de Onda e que depende de inúmeros parâmetros nem sempre controláveis experimentalmente. Apresenta-se, em seguida, um modelo original de simulação dinâmica do comportamento da coluna e que é constituído por um sistema de equações diferenciais não lineares (parâmetros distribuidos). No entanto, se os débitos forem arbitrados como constantes, o sistema passa a ser linear apesar de não possuir solução analítica evidente (p.e. por transformações integrais). A discretização por diferenças finitas permitiu superar estas dificuldades. Existe uma notável concordância entre os valores experimentais e os previstos no modelo.

Simulação de colunas de destilação de elevada eficiência energética

A destilação é o principal processo de separação usado nas indústrias químicas para separação de misturas líquidas, e apesar de ser um processo muito intensivo energeticamente continua hoje em dia a ser o processo de separação preferencial. O largo impacto dos processos de destilação nos custos de investimento e operação tem impulsionado o desenvolvimento de sistemas complexos de destilação, cuja principal finalidade consiste na otimização energética do processo de separação. No âmbito do estudo da utilização de colunas de elevada eficiência energética para separações multicomponente através da simulação de processos, o presente trabalho foca-se na compreensão, análise e avaliação da Coluna de Petlyuk, sistema de colunas com acoplamento térmico integral (FTCDC), inicialmente introduzida na década de 1940, com um reportado potencial para poupanças energéticas até cerca de 30% relativamente ao sistema tradicional. Para a simulação dos sistemas em estudo foi usado o software de simulação de processos HYSYS, cuja aplicação é desenvolvida e analisada no decurso deste trabalho, o qual é dividido em três partes principais: implementação dos processos em HYSYS para simulação rigorosa em modo de estado estacionário, otimização e simulação dinâmica e controlo. A primeira parte diz respeito ao design dos sistemas e âmbito de aplicação do HYSYS nesta etapa. Através do estudo da separação de mistura ternária equimolar de isómeros de butanol, é apresentado um procedimento de implementação do sistema FTCDC e avaliados os resultados de simulação rigorosa fazendo a sua comparação com os resultados obtidos para o sistema tradicional. A segunda parte é dedicada à otimização de ambos os sistemas, FTCDC e tradicional, cujos estudos se baseiam na simulação em estado estacionário. O elevado número de variáveis de design do sistema FTCDC pode conduzir a tedioso e moroso trabalho de simulação. Neste trabalho são aplicadas metodologias de desenho de experiências na otimização do sistema FTCDC, as quais permitiram uma significativa redução o número de simulações a executar, sendo o procedimento de otimização, baseado no procedimento de desenhos de experiências, desenvolvido ao longo desta parte. Os resultados de otimização de ambos os sistemas são avaliados e comparados. A última parte deste trabalho debruça-se sobre o tema da operabilidade dos sistemas, avaliada através da simulação dinâmica dos mesmos e avaliação da estrutura de controlo implementada. A simulação dinâmica de processos, baseada em modelos dinâmicos, permite melhorar o seu design, otimização e operação, cujos resultados traduzirão de uma forma mais realística o desempenho do processo desenhado, uma vez que estes nunca se encontram num verdadeiro estado estacionário. Nesta parte é apresentado o procedimento de implementação da simulação dinâmica de ambos os sistemas, no qual se inclui a definição da configuração de controlo para cada um destes. Através dos resultados de simulação dinâmica de ambos os sistemas é avaliado o seu design e respetiva estrutura de controlo.

Modelagem e simulação da câmara de reação de uma caldeira a carvão pulverizado

O presente trabalho é parte integrante de um projeto multidisciplinar cujo objetivo é construir uma ferramenta computacional capaz de simular a operação de uma Usina Termoelétrica tratando, especificamente, da câmara de reação (local da caldeira onde se desenvolve a queima do combustível). Foi desenvolvido um modelo matemático dinâmico, pseudo-homogêneo, unidimensional e, com propriedades físicas variáveis, com o objetivo de descrever a conversão do calor gerado pela combust em energia térmica, que está baseado em balanços de massa e energia. Este modelo foi implementado no simulador gPROMS e seus parâmetros foram estimados utilizando o módulo gEST deste software (método da máxima verossimilhança), com base em dados de operação da Usina e dados cinéticos de literatura. Para a simulação dinâmica da câmara de reação foi utilizado um método de integração implícita de passos múltiplos com ordem e passo variáveis. A câmara de combustão foi dividida em três zonas, sendo realizados, em cada zona, os balanços de massa e energia, adotando-se a hipótese de mistura perfeita. O modelo desenvolvido prevê o perfil de temperatura ao longo da câmara, as correntes de saída e a composição dos gases de combustão em cada zona; incluindo uma previsão aproximada da concentração dos poluentes formados, principalmente monóxido de nitrogênio e dióxido de enxofre. Os resultados apresentados mostraram que o modelo prediz satisfatoriamente a composição do gás de combustão na saída da câmara de reação e que a divisão da câmara em diferentes zonas de combustão é apropriada para avaliar a geração e o aproveitamento energético na caldeira.

Ferramentas para análise dinâmica e estrutural de um reboque de linha leve

Neste trabalho é desenvolvida uma metodologia de projeto para identificar as regiões críticas da estrutura de um reboque de linha leve sendo tracionado em pavimentos do tipo rodovia de baixa qualidade e estrada secundária de muito baixa qualidade. Para tanto, são levantados alguns dados experimentais da estrutura, necessários para a aproximação e simulação dinâmica de um modelo simplificado. A excitação da base é realizada por atuadores que simulam as oscilações verticais de um perfil de estrada, a qual é definida de acordo com os estudos realizados por Dodds e Robson (1973). Isto permite a determinação de um histórico de carregamentos das regiões da estrutura do chassi sob a ação das molas da suspensão. Em seguida, é gerado um modelo estrutural simplificado do reboque em elementos finitos, chamado de global, no qual são determinadas as regiões sob ação das maiores tensões. Tendo identificada a região mais crítica da estrutura, é criado um modelo local desta parte, onde se pode observar a distribuição de tensões com mais detalhe, permitindo a identificação dos pontos de concentração de tensões. Desta forma, com a aplicação do método de análise global-local é possível a obtenção de resultados detalhados quanto aos esforços da estrutura com um menor custo computacional.

Simulação computacional de medidas estereológicas em estruturas de metal duro (WC-Co)

This work focuses on the creation and applications of a dynamic simulation software in order to study the hard metal structure (WC-Co). The technological ground used to increase the GPU hardware capacity was Geforce 9600 GT along with the PhysX chip created to make games more realistic. The software simulates the three-dimensional carbide structure to the shape of a cubic box where tungsten carbide (WC) are modeled as triangular prisms and truncated triangular prisms. The program was proven effective regarding checking testes, ranging from calculations of parameter measures such as the capacity to increase the number of particles simulated dynamically. It was possible to make an investigation of both the mean parameters and distributions stereological parameters used to characterize the carbide structure through cutting plans. Grounded on the cutting plans concerning the analyzed structures, we have investigated the linear intercepts, the intercepts to the area, and the perimeter section of the intercepted grains as well as the binder phase to the structure by calculating the mean value and distribution of the free path. As literature shows almost consensually that the distribution of the linear intercepts is lognormal, this suggests that the grain distribution is also lognormal. Thus, a routine was developed regarding the program which made possible a more detailed research on this issue. We have observed that it is possible, under certain values for the parameters which define the shape and size of the Prismatic grain to find out the distribution to the linear intercepts that approach the lognormal shape. Regarding a number of developed simulations, we have observed that the distribution curves of the linear and area intercepts as well as the perimeter section are consistent with studies on static computer simulation to these parameters.

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

Simulação de desempenho luminoso para salas de aula em Natal-RN

The city of Natal has a significant daylight availability, although it use isn’t systematically explored in schools architecture. In this context, this research aims to determine procedures for the analysis of the daylight performance in school design in Natal-RN. The method of analysis is divided in Visible Sky Factor (VSF), simulating and analyzing the results. The annual variation of the daylight behavior requires the adoption of dynamic simulation as data procedure. The classrooms were modelled in SketchUp, simulated in Daysim program and the results were assessed by means of spreadsheets in Microsoft Excel. The classrooms dimensions are 7.20mx 7.20m, with windows-to-wall-ratio (WWR) of 20%, 40% and 50%, and with different shading devices, such as standard horizontal overhang, sloped overhang, standard horizontal overhang with side view protection, standard horizontal overhang with a dropped edge, standard horizontal overhang with three horizontal louvers, double standard horizontal overhang, double standard horizontal overhang with three horizontal louvers, plus the use of shelf light in half the models with WWR of 40% and 50%. The data was organized in spreadsheets, with two intervals of UDI: between 300lux and 2000 lux and between 300lux and 3000lux. The simulation was performed with the weather file of 2009 to the city of NatalRN. The graphical outputs are illuminance curves, isolines of UDI among 300lux and 2000 lux and tables with index of occurrences of glare and to an UDI among 300lux 3000lux. The best UDI300-2000lux performance was evidenced to: Phase 1 (models with WWR of 20%), Phase 2 (models with WWR of 40% and 50% with light shelf). The best UDI300-3000lux performance was evidenced to: Phase 1 (models with WWR of 20% and 40% with light shelf) and Phase 2 (models with WWR of 40% and 50% with light shelf). The outputs prove that the daylight quality mainly depends on the shading system efficacy to avoid the glare occurrence, which determines the daylight discomfort. The bioclimatic recommendations of big openings with partial shading (with an opening with direct sunlight) resulted in illuminances level higher than the acceptable upper threshold. The improvement of the shading system percentage (from 73% to 91%) in medium-size of openings (WWR 40% and 50%) reduced or eliminate the glare occurrence without compromising the daylight zone depth (7.20m). The passive zone was determined for classrooms with satisfactory daylight performance, it was calculated the daylight zone depth rule-of-thumb with the ratio between daylight zone depth and the height of the window for different size of openings. The ratio ranged from 1.54 to 2.57 for WWR of 20%, 40% and 50% respectively. There was a reduction or elimination of glare in the passive area with light shelf, or with awning window shading.

Simulação computacional de medidas estereológicas em estruturas de metal duro (WC-Co)

This work focuses on the creation and applications of a dynamic simulation software in order to study the hard metal structure (WC-Co). The technological ground used to increase the GPU hardware capacity was Geforce 9600 GT along with the PhysX chip created to make games more realistic. The software simulates the three-dimensional carbide structure to the shape of a cubic box where tungsten carbide (WC) are modeled as triangular prisms and truncated triangular prisms. The program was proven effective regarding checking testes, ranging from calculations of parameter measures such as the capacity to increase the number of particles simulated dynamically. It was possible to make an investigation of both the mean parameters and distributions stereological parameters used to characterize the carbide structure through cutting plans. Grounded on the cutting plans concerning the analyzed structures, we have investigated the linear intercepts, the intercepts to the area, and the perimeter section of the intercepted grains as well as the binder phase to the structure by calculating the mean value and distribution of the free path. As literature shows almost consensually that the distribution of the linear intercepts is lognormal, this suggests that the grain distribution is also lognormal. Thus, a routine was developed regarding the program which made possible a more detailed research on this issue. We have observed that it is possible, under certain values for the parameters which define the shape and size of the Prismatic grain to find out the distribution to the linear intercepts that approach the lognormal shape. Regarding a number of developed simulations, we have observed that the distribution curves of the linear and area intercepts as well as the perimeter section are consistent with studies on static computer simulation to these parameters.